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connection,  motor  weight  reduces  usable  payload.  Asymptotic  methods,  scattering,  slender  body  theories  and  CFD  modeling  have  provided  valuable 
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provide  scaling  laws  as  well  as  closed-form  expressions  for  the  interference  not  accessible  solely  from  computation  that  can  be  used  for  interpolation 
and  extrapolation  of  CFD  solutions  as  well  as  efficient  testing  and  design  of  new  flight  vehicles.  This  approach  dramatically  simplifies  trajectory 
predictions  in  which  inertial  and  aerodynamic  forces  are  strongly  coupled.  The  aforementioned  modeling  has  been  carried  out  in  the  framework  of 
the  quasi-steady  approximation.  However,  there  are  cases  when  acceleration  of  the  body  C.G.  and/or  the  pitching  angular  velocity  are  not  small  and 
the  unsteady  effects  become  appreciable.  First  steps  in  mathematical  modeling  of  these  effects  have  been  made  with  emphasis  on  the  body  wave 
drag.  Namely,  theoretical  analysis  of  the  wave  drag  generated  by  a  non -lifting  slender  body  of  revolution  during  its  instantaneous  supersonic  start 
was  carried  out.  Typical  phases  of  the  transient  process  were  identified.  For  each  phase,  analytical  solutions  of  the  flow  potential  and  wave  drag  wert 
obtained.  The  transient  process  was  also  simulated  numerically  by  solving  3-D  Euler  equations.  The  theoretical  wave  drag  is  in  excellent  agreemen* 
with  CFD  in  all  phases  of  the  transient  process.  These  results  provide  a  good  launching  pad  for  theoretical  modeling  of  unsteady  body  motions 
including  the  coupling  between  body  dynamics  and  aerodynamics. 
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Executive  summary 

New  insights  have  been  gained  in  atmospheric  and  space  launch  stage  separation  from  our 
studies  of  the  aerodynamic  interference  between  multiple  bodies  in  supersonic  flow.  Quick 
means  of  estimating  and  controlling  repulsion  or  attraction  lift  associated  with  this  interaction  is 
an  important  enabling  technology  to  size  launch  separation  rocket  motors  to  achieve  the  best 
compromise  between  motor  weight  and  safe  staging.  In  this  connection,  motor  weight  reduces 
usable  payload. 

Asymptotic  methods,  scattering,  slender  body  theories  and  CFD  modeling  have  provided 
valuable  systematic  approximations  schemes  that  advantageously  couple  with  modem 
computational  methods.  Our  theoretical  solutions  for  lift  force  interference  between  multiple 
bodies  give  good  agreement  with  numerical  solutions  and  experimental  data.  These  solutions 
shed  light  on  important  scattering  phenomena  not  previously  recognized  as  relevant  to  this 
problem.  The  analyses  allow  us  to  identify  lumped  nondimensional  parameters  and  provide 
scaling  laws  as  well  as  closed-form  expressions  for  the  interference  not  accessible  solely  from 
computation  that  can  be  used  for  interpolation  and  extrapolation  of  CFD  solutions  as  well  as 
efficient  testing  and  design  of  new  flight  vehicles.  This  approach  dramatically  simplifies 
trajectory  predictions  in  which  inertial  and  aerodynamic  forces  are  strongly  coupled. 

The  aforementioned  modeling  has  been  carried  out  in  the  framework  of  the  quasi-steady 
approximation.  However,  there  are  cases  when  acceleration  of  the  body  C.G.  and/or  the  pitching 
angular  velocity  are  not  small  and  the  unsteady  effects  become  appreciable.  First  steps  in 
mathematical  modeling  of  these  effects  have  been  made  with  emphasis  on  the  body  wave  drag. 
Namely,  we  have  conducted  theoretical  analysis  of  the  wave  drag  generated  by  a  non-lifting 
slender  body  of  revolution  during  its  instantaneous  supersonic  start.  The  analysis  is  based  on 
linear  acoustic  theory  and  asymptotic  techniques.  Typical  phases  of  the  transient  process  were 
identified.  For  each  phase,  analytical  expressions  of  the  flow  potential  and  wave  drag  were 
derived.  A  composite  solution  for  the  wave-drag  coefficient  that  is  uniformly  valid  throughout 
the  transient  process  was  obtained.  The  analytical  solutions  were  plugged  into  a  FORTRAN  code 
providing  quick  calculations  of  the  wave  drag  versus  time.  Numerical  examples  were  generated 
for  a  slender  body  of  half-sine  shape  for  different  Mach  numbers  of  instantaneous  body  start.  It 
was  shown  that  the  initial  wave-drag  coefficient  is  significantly  larger  than  its  asymptotic  level 
relevant  to  steady  flight.  The  transient  time  period  increases  and  tends  to  infinity  as  the  Mach 
number  decreases  and  approaches  M  =  1  (transonic  case).  The  transient  process  was  also 
simulated  numerically  by  solving  3-D  Euler  equations  for  Mach=2  start  of  the  aforementioned 
body.  It  was  shown  that  the  theoretical  wave  drag  is  in  excellent  agreement  with  CFD  in  all 
phases  of  the  transient  process.  These  results  provide  a  good  launching  pad  for  theoretical 
modeling  of  unsteady  body  motions  including  the  coupling  between  body  dynamics  and 
aerodynamics. 

Our  experience  in  low-order  modeling  of  store  separation  helped  Dr.  Malmuth  to  develop  a 
physics-based  predictive  tool  for  quick  estimates  of  the  body  trajectory  in  the  case  of  separation 
from  a  rectangular  cavity  to  supersonic  outer  flow.  This  tool  has  been  used  for  parametric  studies 
of  store  trajectories  relevant  to  the  HIFEX  sled  tests.  The  theoretical  predictions  along  with  CFD 
studies  of  Boeing  convinced  the  Air  Force  that  the  HIFEX  store  release  was  not  a  high-risk 
event. 
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1.  Introduction 


Emphasis  on  safe,  high-accuracy,  reliable  parasite  body  separation  will  be  critical  for  parasite 
body  system  upgrades  to  cost-effectively  for  new  aircraft  development  for  the  foreseeable  future. 
Although  the  problem  of  parasite  bodies  release  from  an  aircraft  platform  has  received  much 
attention  as  exemplified  by  Refs.  [1-3],  more  effort  is  needed  to  understand  the  underlying  basic 
physics  and  essential  parameters.  In  a  larger  sense,  this  technology  is  a  subset  of  the  problem  of 
the  interactions  between  moving  bodies  in  all  speed  ranges.  Multi-body  interaction  applications 
include  separation  and  carriage  of  various  stage  vehicles  for  space  missions  and  flight-testing  as 
well  as  crew  escape.  For  this  group,  hypersonic  multistage  vehicle  concepts  of  interest  to  the 
U.S.  Air  Force  frequently  utilize  the  launch  of  a  small  rocket-powered  stage  from  a  large 
subsonic  or  transonic  aircraft  such  as  the  B-52.  Other  examples  are  the  PEGASUS  series  and  the 
Shuttle.  Currently,  multistage  launch  scenarios  are  envisioned  for  future  hypersonic  and  space 
applications.  A  new  thrust  involves  the  use  of  air  breathing  (scramjet  or  turbo-scramjet)  rather 
than  rocket-powered  stages.  These  can  be  associated  with  recoverable  launch  vehicles  such  as 
the  RLV  series  (second  and  third  generation)  being  studied  in  the  Boeing  CRI  program.  Other 
applications  include  the  Boeing  Delta  IV  series. 

The  payload  mass  fraction  can  be  drastically  affected  by  sizing  of  the  separation  motors.  Overly 
conservative  estimates  requiring  extra  thrust  will  give  a  larger  weight  penalty  associated  with 
more  motor  thrust  for  stage  separation.  A  smaller  (optimistic)  separation  motor  design  solution 
will  increase  payload  with  the  risk  of  re-contact  and  reduced  safety.  The  correctness  of  these 
decisions  depends  on  our  understanding  of  the  interaction  of  the  fluid  dynamics  with  the 
dynamics  of  the  motion. 

Rather  than  emphasizing  speed,  efficiency  and  accuracy  of  production  CFD  codes,  we  stress 
computationally  non-intense  PC  application  to  understand  and  harness  (through  new  flow  control 
concepts)  the  important  physical  mechanisms  that  pervade  a  wide  class  of  parasite  body  and 
stage  separation  phenomena.  Another  emphasis  is  to  formulate  representative  unit  problems  and 
identify  new  interesting  mathematical  techniques  to  solve  these  problems  and  combine  with  them 
current  large-scale  numerical  techniques.  On  the  other  hand,  more  realism  is  being  introduced  as 
our  effort  progresses. 

Besides  the  research  thrust  of  stage  separation  embodied  in  the  effort  discussed  herein,  another  is 
a  tie-in  with  current  efforts  regarding  optimizing  parasite  bodies-bay  acoustic  mitigation 
measures  with  stage  trajectory  characteristics,  (see  [4]).  Current  control  measures  include 
alternate  jet  blowing  and  actuators  to  reduce  the  amplitudes  of  the  noise  power  from  parasite 
bodies-bay  cavities.  Flight  and  ground  tests  are  underway  to  investigate  the  effects  on  parasite 
body  trajectories.  The  physics  and  time  scales  of  the  mitigation  measures  and  their  possible 
coupling  with  the  dynamic  characteristics  of  the  separation  bodies  needs  to  be  better  understood. 
Our  models  provide  useful  insight  into  this  interaction.  Examples  of  our  work  are  in  Refs.  [5-16]. 

2.  Major  accomplishments 

A  major  thrust  in  FY  2005-2006  was  a  mathematical  modeling  of  lift  and  pitching  moment 
associated  with  the  two-body  interaction.  The  general  mathematical  ideas  were  discussed  in  [16] 
and  key  unit  problems  are  shown  in  Figure  1.  Special  emphasis  was  on  Problems  1,  2  and  4 
relevant  to  scattering  of  the  parent-body-induced  shock  wave  by  the  parasite  body. 
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In  FY  2005  we  studied  the  interaction  between  the  wedge-induced  shock  and  a  cylinder  using  the 
theory  of  scattering  (Problem  1)  [15].  The  lift  force  associated  with  multi-scattering  of  the  wing- 
induced  shock  wave  by  a  slender  body  of  revolution  was  analyzed  using  linearized  supersonic 
theory,  scattering  theory  [17,18]  and  asymptotic  methods  [19-21].  Although  acoustical  and 
electromagnetic  analyses  such  as  those  just  cited  are  well  known,  their  application  to  supersonic 
shock-wave  scattering  problems  has  not  received  significant  attention.The  local  and  integral  lift 
coefficients  were  obtained  in  simple  analytical  forms  convenient  for  quick  calculations  of  the 
aerodynamic  loads.  These  solutions  can  be  applied  to  the  moderate  supersonic  (Mach  numbers 
from  1.2  to  3)  multi-body  interaction  problem  for  crosscheck  with  other  computational  or 
engineering  methods. 


Problem  1:  Plane 
shock-parasite  body 
interaction 


Problem  2: 
parasite  body- 
induced 
wave  train 


Problem  3:  parasite  body-induced 
wave-train  including  expansion  waves 


Problem  5:  Wave  train  between  two 
wedges 

supersonic  stage  separation. 


The  local  analytical  solution,  which  provides  detailed  distribution  of  the  lift-force  coefficient  on 
the  length  scale  ~  body  radius,  is  shown  by  the  black  line  in  Figure  2.  The  global  solution  ignores 
local  details  and  gives  the  step-function  distribution  shown  by  the  red  line.  It  was  found  that  the 
total  lift  is  predominantly  generated  by  the  first  scattering.  We  conduced  series  of  Euler  CFD 
calculations  for  Mach=2  free  stream  (Figure  3),  and  showed  that  the  total  lift-force  coefficient 
predicted  by  CFD  (red  symbols)  agrees  well  with  the  theoretical  solution  (black  line).  The 
analytical  solution,  which  is  inaccessible  from  purely  numerical  methods,  shows  the  excellent 
complementarity  of  the  numerics  and  the  pen-and-paper  methods.  On  the  other,  the  numerics 
shows  accessibility  to  nonlinear  effects  not  modeled  by  the  linear  theory. 

We  also  showed  that  the  analytical  solution  agrees  satisfactory  with  the  experimental  data  [22] 
(Figures  4  and  5)  obtained  for  the  blunt-base  parabolic  body  of  revolution  in  the  flow-field 
induced  by  a  circular-arc  wing  in  the  Mach=2  free  stream,  despite  the  fact  that  the  body  radius 
varies  in  the  interaction  region  and  the  wing  thickness  is  not  very  small  (the  wing  thickness  ratio 
is  1/6).  This  adds  more  confidence  in  practicability  of  the  theoretical  model  [15]. 
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P  ft  X 


Global  solution  cl(x)=4sH(x-x0) 


Figure  2  Lift  force  CL  induced  by  scattering  of  the  parent-body-induced  shock  wave  by  the 
parasite  body. 


Freestream  Mach=2 
Body  nose  length  l=20a 
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Body  is  completely  behind  shock 
(uniform  upwash) 
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wedQe  '  '  base  '  nose'1  nose 


0.16 


0.14 


—  theory 
9  CFD 
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0.00 


-15  -10  -S  0  5/  10  15  20  25  30  35 
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CL=4ct 


Body  nose  is  completely 
upstream  from  shock 
(shock-cylinder  interaction) 


x  =  distance  from  nose  tip  to 
point  where  shock  crosses 

body  axis _ 


Shock  crosses  body  nose 


Figure  3  The  total  lift  force  predicted  by  theory  (black  line)  agrees  well  with  the  Euler  CFD 
solution  (symbols  and  panels  showing  the  pressure  field). 


6 


Ron  View 


Wing  Section 


Figure  4  Schematic  layout  of  test  setup  and  models  (scanned  from  Figure  1  of  Ref.  [22]). 


Figure  5  Comparison  of  theoretical  (solid  line)  and  experimental  (symbols)  distributions  of  the 
normal  force  induced  by  the  wing  shock  impinging  on  the  body  at  x  ps  0.4 . 

In  FY  2006,  we  treated  the  case  of  two-body  interference  schematically  shown  in  Figure  6  [23]. 
For  this  configuration,  the  cross  sectional  shape  of  both  bodies  is  a  function  of  the  axial 
coordinate,  and  bodies  generate  3-D  flow  disturbances  with  appreciable  streamwise  variations 
over  the  body  length. 
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3-D  problem:  /?6^-(s^+a^)= o 

dn<j>-Q  on  parent-body  surface 

dn(f>=0  on  parasite-body  surface 

Figure  6  Problem  of  the  two-body  interaction. 


The  problem  is  analyzed  using  two  essentially  different  approaches: 

•  Approach  1  provides  a  global  solution.  The  problem  is  solved  using  asymptotic  methods  as 
well  as  the  basic  idea  of  slender  body  theory  -  replace  bodies  by  distributions  of  sources  and 
doublets 

•  Approach  2  is  based  on  the  shock-wave  scattering  theory.  In  this  connection,  the  theoretical 
model  [15]  developed  in  2005  was  extended  to  the  case  of  two  bodies  of  revolution 

These  two  approaches  complement  each  other:  the  first  provides  a  global  structure  of  the  lift 
force  distribution,  the  second  provides  flow  details  in  local  regions  associated  with  rapid  changes 
of  incoming  (parent-body-induced)  flow  field  in  the  vicinity  of  the  parasite  body  such  as  shock- 
body  interaction  regions. 

We  started  with  Approach  1  and  considered  a  parasite  body  in  the  parent-body-induced 
nonuniform  supersonic  flow  (Figure  7).  For  large  vertical  distance  between  bodies  (h/ a  ~  e~l), 
the  parent-body-induced  disturbance  is  predominantly  two-dimensional  (to  relative  accuracy 
0(s2))  in  the  vicinity  of  parasite  body.  This  potential  gives  the  vertical  velocity  on  the  parasite 
body  axis  (see  the  first  box).  The  vertical  velocity  perturbation  is  treated  as  the  induced  angle  of 
attack.  The  parasite  body  is  replaced  by  line-distributions  of  sources  and  doublets.  In  the  vicinity 
of  parasite  body  axis,  the  radial  velocity  (induced  by  parasite  body)  has  the  form  shown  in  the 
second  box.  From  boundary  conditions  on  the  parasite  body  surface  we  obtain  intensities  of 
sources  and  doublets  (see  the  third  box).  Then  we  calculate  the  pressure  coefficient  on  the 
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parasite  body  surface  and  extract  its  asymmetric  part,  which  contributes  to  the  lift  force  (see  the 
fourth  box).  Integrating  the  asymmetric  part  of  pressure  coefficient  over  the  parasite  body 
surface,  we  obtain  the  axial  distribution  of  the  lift-force  coefficient  (see  the  pink  box). 


Figure  7  Major  steps  of  Approach  1  providing  the  global  solution. 


To  illustrate  basic  features  of  the  two-body  interaction  we  analyzed  the  interaction  between  two 
identical  bodies,  which  have  a  sine-shaped  nose  of  length  ln  and  a  long  cylindrical  afterbody. 

The  vertical  distance  between  bodies  was  fixed,  h  =  ln,  whereas  the  streamwise  distance  between 
the  nose  tips  was  variable.  Figure  8  shows  the  normalized  lift-force  coefficient  as  a  function  of 
x0  -  streamwise  distance  from  the  parasite  body  tip  to  the  point  where  the  parent-body-induced 
shock  crosses  the  parasite-body  axis. 

The  yellow  frame  in  Figure  8  indicates  the  case,  when  the  parent-induced  disturbance  interacts 
with  the  cylindrical  afterbody  only.  This  case  has  been  studied  in  details  using  Approach  2  based 
on  the  scattering  theory.  Basic  steps  of  this  analysis  are  shown  in  Figure  9.  The  analysis  is  started 
with  replacing  of  the  parent  body  by  the  line-source  distribution.  The  Laplace  transform  of 
corresponding  perturbation  potential  contains  modified  Bessel  function  K0  (see  the  first  box). 

Asymptotic  expansion  of  this  transform  for  large  argument  gives  a  solution,  which  can  be  treated 
as  a  plane  wave  radiating  the  parasite-body  cylindrical  surface.  The  interaction  of  each  plain 
wave  with  the  cylinder  is  solved  analytically  using  methods  of  the  scattering  theory.  Solution  of 
this  unit  problem  provides  the  pressure  distribution  on  the  parasite  body.  Integrating  this  pressure 
over  the  parasite  body  surface  we  obtain  Fourier  components  of  the  local  lift  coefficient  from 
each  elementary  wave  (see  the  second  box).  The  inverse  Laplace  transform  gives  the  local  lift 
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coefficient  per  unit  length  (see  the  third  box).  Finally,  we  obtain  the  axial  distribution  of  lift- 
force  coefficient  shown  in  the  pink  box. 


parant  body 

0<x0<lnose:  shock  crosses  body  nose 
Negative  CL  indicates  attraction 


x0>lnose:  shock  crosses 
cylindrical  afterbody 
CL=0,  no  lift 


i  i 

=  2  ja(x)S'(x)dx  +  4  ja’(x)S(x)dx  = 

I 

=  4a(l)S(l)-2ja(x)S,(x)dx 

Gl  s  CL  tgtaI !  e2  -  normalized  lift-force  cocfftcieni 


Mach=2,  h=l  =1 

2  - - — 


xo<0:  shock  upstream  from  nosetip 
Positive  CL  indicates  repulsion 


Figure  8  Lift  due  to  the  interaction  between  two  identical  bodies  of  revolution  (Approach  1). 
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Figure  9  Basic  steps  in  solving  the  scattering  problem  for  the  shock  scattering  on  the  cylindrical 
afterbody  (Approach  2). 
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To  verify  aforementioned  theoretical  models,  we  carried  out  3-D  Euler  simulation  of  this 
problem  (Figure  10).  For  two  identical  bodies  with  aligned  nosetips,  the  problem  is  identical  to 
that  for  one  body  over  a  horizontal  plane  at  the  vertical  distance  hi 2  from  the  body  axis.  The 
computational  domain  for  this  configuration  is  shown  in  the  left  bottom  comer.  The  pressure, 
temperature  and  Mach  number  distributions  in  the  symmetry  plane  as  well  as  the  pressure 
distribution  on  the  body  surface  indicate  a  complex  structure  of  induced  and  reflected 
disturbances. 


Ptrtnl  tody 

xo>,nose:  shock  crosses 
cylindrical  afterbody 

CL=0,  no  lift 


Computational  domain 


iiia 

jg| 

ISii 

m 


plain  wall 


plain  wall 


Pressure  field  Temperature  field 


Mach  number  field 


Mach=2,  nose  length  lnose=18a,  h=16a 


Figure  10  3-D  Euler  CFD  solution  for  two  identical  slender  bodies  of  revolution. 


Figure  11  compares  the  theoretical  distributions  of  CL(x )  with  the  CFD  solution.  The  blue  line 

shows  the  global  solution  obtained  with  Approach  1  (see  Figure  7).  The  black  line  corresponds  to 
the  local  solution  obtained  with  Approach  2  based  on  the  scattering  problem  (see  Figure  9).  The 
red  symbols  show  the  Euler  CFD  data.  Despite  3-D  complexity  of  two-body  interference,  the 
leading-order  analytical  solutions  agree  with  the  CFD  to  0{s )  . 

We  also  considered  the  scaling  issue  of  the  two-body  interaction  (Figure  12).  The  analytical 
forms  of  global  and  local  asymptotic  solutions  allowed  us  to  formulate  the  scaling  law  (see  the 
pink  box).  Robustness  of  this  scaling  was  estimated  using  parametric  calculations  of  the  lift-force 
coefficient  as  a  function  of  x0  (streamwise  distance  from  the  nose  tip  to  the  point  where  the 

shock  crosses  the  body  axis)  at  various  vertical  distances  between  two  identical  bodies  (see  the 
left  plot)  and  various  Mach  numbers  (the  right  plot).  It  is  seen  that  these  dependencies  quickly 
collapse  to  a  universal  curve  as  the  separation  distance  h  and/or  the  freestream  Mach  number 
increases.  Using  this  scaling  it  is  feasible  to  interpolate/extrapolate  CFD  solutions  and 
dramatically  reduce  the  computational  effort  that  is  of  critical  importance  for  quick  assessments 
of  aerodynamic  coefficients  at  various  time  moments  during  the  separation  process. 
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Mach=2,  nose  length  lnose=18a,  h=16a 

Figure  11  Comparison  of  theoretical  and  CFD  solution  for  the  lift  force  coefficient  distribution 
along  the  parasite  body  axis. 
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Figure  12  Scaling  issue  for  interference  between  two  identical  bodies  of  revolution. 
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The  aforementioned  modeling  has  been  carried  out  in  the  framework  of  the  quasi-steady 
approximation.  However,  there  are  cases  when  acceleration  of  the  body  center  of  gravity  and/or 
the  pitching  angular  velocity  are  not  small  and  the  unsteady  effects  become  appreciable.  In  FY 
2007,  we  started  a  mathematical  modeling  of  these  effects.  As  a  first  step  we  considered  the 
wave  drag  associated  with  unsteady  motion  of  a  slender  body  of  revolution  as  shown  in  Figure 
13. 


Figure  13  Slender  body  of  revolution  moving  along  x  -axis  with  speed  U ( t ) . 

It  is  assumed  that  the  body-induced  disturbance  is  small  compared  with  speed  of  sound,  and  the 
analysis  is  based  on  the  linearized  3-D  Euler  equations.  The  disturbance  velocity  potential 
$(x,y,z,t)  is  governed  by  the  wave  equation 


vV 


(i) 


where  ax  is  speed  of  sound  in  undisturbed  fluid.  A  slender  body  moving  along  the  x  -axis 
induces  unsteady  acoustic  sources  whose  intensity  Q(x,t )  is  coupled  with  the  body  shape  and 
motion  as  Q(x,  t)  =  dAh  (x,  t)/  dt ,  where  Ab(x,t)  =  7rrb(x,t)  is  the  body  cross-sectional  area.  For 
the  body  of  constant  shape  and  length  L  moving  with  the  velocity  U (t) ,  we  have 

t 

Ab(x,t )  =  Ab(X) ,  where  X  =  x+  ,  0  <  X  <  L ,  the  solution  is  expressed  in  terms  of  the 

0 

retarded  potential  as  [21,24] 


,,  1  %  U(t)4(tj)  ,e 

</>{x,y,z,t)  =  -—  J  .  =d^, 

r  =  ^  +r  - ,  tj  =  ^+  j[/(r')Jr',  r  =  yfy^  +z2  . 


(2) 

(3) 


The  pressure  distribution  on  the  body  surface  is  calculated  using  the  near-field  potential 
<j)(x,r  — »  0,t)  and  the  unsteady  Bernoulli  equation  in  the  slender-body  approximation. 

We  considered  an  instantaneous  acceleration  of  the  body  to  a  certain  constant  speed  U  at  the 
initial  time  t  =  0  and  analyzed  the  wave  drag  D(t)  during  the  transient  process.  Figure  14  shows 
trajectories  (black  lines)  of  the  body  leading  and  trailing  edges,  the  typical  regions  A,  B  and  C 
associated  with  phases  of  the  transient  process  as  well  as  characteristics  (red  dashed  lines)  and 
end  points  x12  typifying  each  phase.  Analytical  expressions  for  the  wave  drag  D(t)  were 

derived  for  each  of  the  aforementioned  phases.  Using  these  expressions  D(t)  was  calculated  and 
compared  with  unsteady  Euler  CFD  solution  for  the  body  of  shape  rh(X)  =  rmax  sin(^Y /  L) , 
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rmax  =  0.0 1Z,  at  the  Mach  number  M  =  U/aa0=2  (Figure  15).  The  theoretical  solution  is  in 

excellent  agreement  with  CFD.  This  analysis  provides  a  launching  pad  for  theoretical  modeling 
of  the  unsteady  body  motion  accounting  for  the  coupling  between  body  dynamics  and 
aerodynamics. 

Details  of  this  study  are  given  in  Appendix  A. 


Figure  14  Characteristics  lines  and  boundaries  for  a  slender  body  accelerating  to  a  supersonic 
speed  U  at  the  initial  time  moment  t  =  0 ;  here  a  =  ax ,  at  t  =  0  the  body  leading 

edge  xle  (0)  =  0  and  the  trailing  edge  xte  (0)  =  L . 


T=ta/L 


T=ta/L 


Figure  15  Comparison  of  theoretical  and  CFD  history  of  drag  coefficient  CD(T), 
Cd(T)=D(T)/(a^P,U‘/2),  T  =  taJL. 
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Our  experience  in  low-order  modeling  of  store  separation  helped  Dr.  Malmuth  to  develop  a 
physics-based  model  for  quick  estimates  of  the  body  trajectory  in  the  case  of  separation  from  a 
rectangular  cavity  to  supersonic  outer  flow.  The  cavity  is  assumed  to  be  deep  with  the  shear  layer 
bridging  the  front  and  rare  cavity  lips  (open  cavity  flow).  3 -DOF  equations  for  the  body  motion 
includes  aerodynamic  loads  (drag,  lift  and  pitch  moment),  which  are  evaluated  analytically  and 
numerically  for  the  three  pahses  of  separation: 

1 .  Body  moves  incide  cavity 

2.  Body  crosses  shear  layer 

3.  Body  moves  outside  cavity  in  supersonic  flow 

In  Phase  1,  the  body  experiences  aerodynamic  loads  due  to  the  flow  circulation  inside  the  cavity 
and  due  to  the  body  motion.  Estimates  based  on  viscous  flow  physics  showed  that  these  loads  are 
of  ther  order  of  0(Re~1)  +  0(p2) ,  where  ReL  is  Reynolds  number  based  on  free-stream 
parameters  and  L  =forebody  length+cavity  length,  (3  =  ve/U  is  ratio  of  the  body  vertical 
velocity  to  the  ffeestream  speed.  For  many  practical  cases,  these  parameters  are  small,  and  the 
aerodynamic  loads  incide  cavity  can  be  ignored.  The  body  dynamic  equations  have  a  simple 
analytical  solution  relevant  to  the  free  drop  in  vacuum. 

In  Phase  2,  the  shear  layer  is  assumed  to  be  much  thiner  than  the  body  radius  and  it  is 
approximated  by  a  straight  line  bridging  the  front  and  rare  cavity  lips.  In  addition,  the  body 
thickness  ratio,  £ ,  is  treated  as  a  small  parameter  (slender  body  approximation).  With  these 
assumptions,  the  aerodynamic  loads  are  expressed  in  analytical  forms  using  the  method  of 
matched  asymptotic  expansions. 

In  Phase  3,  the  equations  for  the  body  trajectory  are  solved  analytically  in  the  first-oreder 
approximation  with  respect  to  small  parameters  £  and  ft .  This  solution  captures  an  osccilatory 
behavior  of  the  pitch  angle  that  depends  on  the  aerodynamic  coefficient  Cma  (the  pitch  moment 
coefficient  derivative  with  respect  to  angle  of  attack). 


Figure  16  Trajectories  at  various  Cma ,  free-stream  Mach  number=2. 
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As  an  example,  Figure  16  shows  the  body  trajectories  at  different  values  of  Cma .  Flow  and  body 
parameters  are  relevant  to  the  HIFEX  sled  test  [25],  As  Cma  tends  to  zero,  the  pitch  angle  0 
evolves  from  oscillatory  to  monotonic  behavior  versus  x .  Ultimately,  0  becomes  a  linear 
function  of  time  at  Cma  =  0  that  gives  monotonic  decreasing  of  0  with  x  (the  green  line).  Note 

that  analytical  solutions  (dashed  lines)  are  close  to  numerical  ones  (solid  lines)  in  all  cases 
considered. 

Similar  parametric  calculations  showed  that  models  would  hit  the  diverter  and  would  not  damage 
the  test  track  (see  Figure  17).  This  analysis  along  with  other  numerical  studies  of  Boeing 
convinced  the  Air  Force  that  the  FIIFEX  store  release  was  not  a  high-risk  event  and  gave  Boeing 
the  go-ahead  for  the  sled  mission.  Details  of  these  studies  can  be  given  by  Bill  Bower  (Boeing). 


Mach  2  rocket  sled  configuration' 
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Figure  17  Physics-based  model  helps  to  estimate  store  trajectories  and  shows  that  models 
released  from  the  HIFEX  bay  do  not  damage  the  test  track  (HIFEX  program). 


3.  Summary  and  impact  of  effort 

In  summary,  new  deductive  asymptotic  framework  has  been  developed  to  treat  the  interaction 
between  two  supersonic  bodies.  Compact,  closed-form  analytical  solutions  were  obtained  for  the 
two-body  supersonic  interference.  These  solutions  are  ideally  suited  for  preliminary  design, 
scaling  and  separation  motor  payload  safety  compromise  that  are  not  available  from  purely 
numerical  solutions.  They  also  provide  a  launching  pad  for  systematic  refinements  and  coupled 


16 


dynamic/aerodynamic  modeling  needed  for  staging  tradeoff  studies.  Our  theoretical  model 
indicated  that  there  are  two  physical  mechanisms  that  control  the  interference  flow  field: 

•  Streamwise  stratification  of  the  transverse  velocity  imposed  from  one  body  on  the  other 
(mutual  induction) 

•  Shock  scattering  from  one  body  incident  on  the  other 

The  aforementioned  analysis  provides  for  the  first  time  an  approximate  but  coherent  picture  of 
the  wave  train  structure  occurring  in  multi-body  interference  problems  such  as  supersonic  stage 
separation  and  flight  crew  escape.  This  is  particularly  important  for  internal  carriage  and  cocoon 
release  where  shear  layers  bounding  cavities  exist.  Moreover,  this  work  provides  a  launching 
pad  for  the  effort  to  obtain  more  understanding  of  the  near  field  as  well  as  increased  prediction 
accuracy.  This  will  strongly  impact  ejection  system  design  and  safe  separation  flow  control. 

Bill  Bower  of  Boeing  is  using  our  store  prediction  and  control  expertise  for  HIFEX  and  SEAR 
programs. 

In  the  case  of  superlight  bodies,  when  acceleration  of  the  body  C.G.  and/or  the  pitching  angular 
velocity  are  not  small,  the  unsteady  effects  become  appreciable.  First  steps  in  mathematical 
modeling  of  these  effects  have  been  made  with  emphasis  on  the  body  wave  drag.  Namely, 
theoretical  analysis  of  the  wave  drag  generated  by  a  non-lifting  slender  body  of  revolution  during 
its  instantaneous  supersonic  start  was  carried  out.  Typical  phases  of  the  transient  process  were 
identified.  For  each  phase,  analytical  solutions  of  the  flow  potential  and  wave  drag  were 
obtained.  The  transient  process  was  also  simulated  numerically  by  solving  3-D  Euler  equations. 
The  theoretical  wave  drag  is  in  excellent  agreement  with  CFD  in  all  phases  of  the  transient 
process.  These  results  provide  a  solid  foundation  for  theoretical  modeling  of  unsteady  body 
motions  including  the  coupling  between  body  dynamics  and  aerodynamics. 
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Appendix  A:  Wave  Drag  of  a  Slender  Non-Lifting  Body  during 
Transient  Process  of  Instanteneous  Supersonic  Start 


1.  Acoustic  theory  for  a  slender  non-lifting  body 

Consider  a  slender  body  of  revolution  moving  with  an  unsteady  speed  U(t)  at  zero  angle  of 
attack  along  the  x  -axis  as  shown  in  Figure  A1 . 


Figure  A1  Slender  body  of  revolution  moving  along  the  x  -axis  with  speed  U (t) . 

Let  q (x,y,z,t)  -  velocity  disturbance  induced  by  the  body  motion  in  gas  at  rest.  It  is  assumed 
that  this  disturbance  is  small:  |q|  /  ar  « 1 ,  where  ax  is  speed  of  sound  in  undisturbed  gas  at  the 

infinity.  The  body-induced  disturbances  are  treated  as  acoustic  disturbances  and  governed  by  the 
linearized  3-D  Euler  equations  of  acoustic  theory 


^  +  Vpq  =  0, 

ct 
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f+"-vq. 
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pr  pi 


=  -vp. 


(l.la) 

(1.1b) 

(1.1c) 


where  y  is  specific  heat  ratio.  For  the  disturbance  velocity  potential  <f>(x,y,z,t) : 
q(x,  y,  z,  t)  =  V </) ,  we  obtain  the  wave  equation 


VV-44 t  =  o,ai=r—- 

tfoc  St  A, 

Disturbances  of  pressure,  p  =  (P-  Ps )  /  Px ,  and  density,  J  =  (p-  )  /  prf ,  are  expressed  as 

p=rj=-&L?£. 

P„  dt 
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A  pointed  spherical  source,  which  is  located  at  (x,y,z)  =  (0,0,0) ,  generates  a  spherical  wave 
with  the  potential 
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This  spherical  wave  induces  the  radial  velocity 
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(1.5) 


The  first  term  of  (1.5)  gives  the  near  field  and  the  second  -  the  far  field.  For  R—> 0 ,  the  mass 
flux  perturbation  (in  units  of  px )  is 


lim4 xR2qR(R,t)  =  -Anf{f)  =  Q(t) , 

x — >oo 

where  Q(t)  -  the  source  intensity.  For  the  impulse  source,  Q(f)  =  S(t) ,  we  obtain 
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(1.7)  is  solution  of  the  problem 
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The  body  induces  unsteady  acoustic  sources  with  the  source  intensity  depending  on  the  body 
shape  and  motion.  If  Q{pc,t)  is  the  source  intensity  at  a  point  (x,t) ,  then  in  accord  with  (1.7) 
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r  =  -y/y2  +z2  . 


(1.10b) 

The  flow  potential  at  ( x,r,t )  is  induced  by  sources  that  were  active  at  early  time  instants 

between 
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z  =  t - .  These  sources  are  located  on  the  hyperbola  z  =  t- 

@CD  "OO 

trajectories  of  the  leading  and  trailing  edges  of  the  body  as  schematically  shown  in  Figure  A2. 


li 


Figure  A2  Influence  domain  for  accelerating  body. 


Following  Ref.  [1]  we  consider  asymptotic  behavior  of  the  integral  (1.10a)  for  r—>  0.  In  this 
limit,  the  dominant  contribution  is  associated  with  the  vicinity  of  £  =  x .  The  integral  (1.10a)  is 
expressed  as  a  sum  of  three  integrals 

0(x,r,t)  =  <t>I+0!I+0I!I,  (1-11) 
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With  the  substitution  £  =  x + r  sinh  a ,  the  second  integral  is  written  as 
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Consider  a  distinguished  limit  s(r )  — »  0 :  r /£—>  0  for  r —>  0 .  Then  (1 . 15)  is  expanded  as 
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Consider  the  first  integral  (1.12).  Integrating  by  parts  we  get 
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where  xx(x,t)  =  xle(x,0,t)  is  shown  in  Figure  A2  (see  intersection  of  the  dashed  line  with  the 
leading-edge  trajectory). 

In  a  similar  way 
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where  x,  (x,/)  =  x(l,(x,0,0  is  intersection  of  the  dashed  line  with  the  trailing-edge  trajectory 
(Figure  A2).  The  final  expression  for  the  flow  potential  in  the  vicinity  of  slender  body  is 


IV 


-47r<f>(x,r,t)  =  2Q(x,t)  log— 

r 

+0  log(x  -xl)  +  Q2  log(x2  -  x) 


+  Jlog (x-^)|a+^-fi|  d^-  j log(£ - x) ja- ~ Q I  ^ 


(1.19) 


+0(r2  logr) 


where 


0i=0 


x-x, 


\ 


a, 


.  02=0 


{/},=/ 


#.<-1 


oc  / 

X 


X2,t- 


x2  —  X 


ao  y 


it 


a 


oo  y 


(1.20) 

(1.21) 

(1.22) 


Xj  (x,  t )  =  xte  (x,  0,  t )  ,  x2  (x,  0  =  xte  (x,  0,  t )  . 

Note  that  the  solution  (1.19)  coincides  with  the  solution  (2)  of  Ref.  [2], 

Now  we  can  couple  the  source  intensity  Q(x,  t)  with  the  body  motion  by  imposing  the  boundary 
condition  on  the  body  surface.  For  the  body  shape 

B(x,r,t)  =  0  =  r-rb(x,t),  (1.23) 

the  boundary  condition  reads 
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where  rb  is  body  radius.  In  the  first-order  approximation,  we  obtain 

.  ,  .  dr. 

<Pr(x,rb,t)  =  - K 
dt 

Using  (1.19),  which  was  derived  for  small  r ,  we  get 

dAb(x,t) 
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dt 


(1.24) 
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where  Ab(x,t)  -  7irb(x,t)  is  the  body  cross-sectional  area. 

For  a  body  of  constant  shape  and  length  L  moving  with  velocity  U (t) ,  we  have 

Ab(x,t)  =  Ab(X),0<X<L,  (1.28a) 
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In  this  case 
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Substitution  of  (1.29)  into  (1.10a)  gives 

1  "J- 

<f>(x,y,z,i)=  f 
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Note  that  this  expression  coincides  with  (1 .3. 1 1)  of  [3]. 


2.  Analysis  of  the  wave  drag 


The  wave  drag  can  be  also  calculated  using  the  near-filed  potential  <f>(x,r  — >0,t)  given  by 
(1.19),  and  the  unsteady  Bernoulli  equation.  For  a  slender  body  (r  — >  0 ),  the  Bernoulli  equation 
reads 


P  =  P~P* 


-Ac 


dt  2  v  dr  ) 


(2.2) 


and  the  wave  drag  is  calculated  as 

D{t)  =  ]pw(X,t)A'b(X)dX ,  (2.3) 

0 

where  X  is  measured  from  the  body  leading  edge,  and  the  pressure  pw  is  calculated  on  the  body 
surface. 

Herein  we  consider  an  instantaneous  acceleration  of  the  body  to  a  certain  constant  speed  U  at 
the  initial  time  moment  t  =  0  and  analyze  the  wave  drag  D(t)  during  the  transient  process. 

Hereafter  we  consider  the  supersonic  case  U  >  ax  following  the  plan: 

1.  Identify  typical  phases  of  the  transient  process  and  determine  the  end  points  x12  (see  Figure 
A2)  associated  with  the  basic  expression  (1.19)  for  the  flow  potential 

2.  Using  (1.19)  calculate  the  near- filed  flow  potential  for  each  phase 

3.  Using  (2.2)  calculate  the  wall  pressure  pw{X,t) 

4.  Using  (2.3)  calculate  the  wave  drag  D(t ) 


vi 


2.1  Structure  of  the  influence  domain 

Figure  A3  shows  trajectories  of  the  leading  and  trailing  edges,  the  typical  regions  A,  B  and  C 
associated  with  phases  of  the  transient  process  as  well  as  locations  of  characteristics  and  end 
points  xl  2  for  each  phase. 

1.  If  the  point  (x,t)  lies  in  Region  A,  the  end  points  x12  are  at  the  body  leading  and  trailing 
edges.  The  wall  pressure  pw(x,t)  =  pwA(X)  should  correspond  to  a  steady  flight  with 
supersonic  speed  U . 

2.  If  the  point  (x,  t)  lies  in  Region  B,  then  Xj  is  at  the  leading  edge,  while  x2  lies  on  the  body 
axis  at  t  =  0 .  Here  the  wall  pressure  pw(x,t )  =  pwB{X,t )  is  unsteady. 

3.  In  Region  C,  the  end  points  x12  are  on  the  body  axis  at  t  =  0,  and  the  wall  pressure 
Pjx,t)  =  Pw,c(X’t)  is  unsteady. 


x1  and  x2  on  leading  edge  x1  on  leading  edge,  x2  -  at  f=0  x1  and  x2  at  t=Q 

Figure  A3  Characteristic  lines  and  boundaries  for  a  slander  body  accelerating  to  a  supersonic 
speed  U  at  the  initial  time  moment  t  =  0 ;  here  a  =  aaz,  at  t  =  0  the  body  leading 

edge  xk  (0)  =  0  and  the  trailing  edge  xte  (0)  =  L . 


For  t  >  L  /(U-ax) ,  the  region  A  covers  the  total  body  length;  i.e.,  the  wall  pressure  corresponds 

to  the  steady  flight  allover  the  body  length.  The  transient  process  is  over  and  the  wave  drag 
should  be  given  by  the  classic  formula  of  von  Karman  and  Moore  [4] 
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Finally  we  get  for  the  near-filed  potential  at  (x,  t)e  A 

X 

2k 


A"(cr)d<j  +  0(r 2  log  r ) 


<f>A(x,r,t)  =  -^-  Jlog ^  j = 


(2.14) 


Note  that  this  potential  coincides  with  the  well  known  near-field  potential  for  steady  supersonic 
flight  (see,  for  example,  Ref.  [5]). 


2.2.2  Potential  in  Region  B 
In  this  region  the  end  points  are 
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Since  A'h  (0)  =  0  for  body  with  a  sharp  nose,  we  have 

Ql=0,Q2=UA'h(x  +  aJ).  (2.19) 

Then,  (1.19)  gives 
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Equation  (2.20)  can  be  written  as 
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Using  the  new  variable  X2 ,  we  get 
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Finally  we  obtain  for  (x,t)  e  B 
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2.2.3  Potential  in  Region  C 


Coordinates  of  the  end  points  are 

*1C  =X~aJ,  X2C=X  +  aJ- 
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The  sources  are 
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Then,  (1.19)  gives 
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This  expression  can  be  written  as 
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Using  the  new  variables  Xx  and  X2 ,  we  obtain 
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Then  the  potential  (2.27)  is  expressed  as 
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Finally  we  get  for  (x,  f)  e  C 


An 


0c(x,r,t)  =  0B(x,r,t)-^-A'b(x-aj)log[(U  +  0/]  +  ^-  j  log (X - cr)  Ab(a)d(r .  (2.28) 
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2.3  Wave  drag  for  t  >  L /(U  -ax) 


Figure  A3  shows  that  for  t  >  L  /(U  -am)  the  body  is  totally  in  Region  A.  Then 
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The  wall  pressure  in  Region  A  is 
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Substituting  (2.32)  into  (2.3)  we  get  the  wave  drag 
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Since  A'b  =  2 Trr.r! ,  we  can  combine  the  first  and  third  terms  as 
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where  we  utilized  the  condition  ^(0)  =  ^  (Z)  =  0 .  The  second  term  is  integrated  by  parts  as 
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Finally  we  obtain  the  classic  formula  of  Karman  and  Moore 
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2.4  Wave  drag  for  L !{U  +  aao)<t<L !{U - ax) 


In  this  case,  Region  A  corresponds  to  0  <  X  <  LA ,  and  Region  B  corresponds  to  LA<  X  <  L , 
where  LA=(U  - a^)t .  The  wave  drag  is 

D{t)=  \pwA(X,t)A'b{X)dX+  J  pwB  (X  ,t)A’h(X)dX  =  DA+  )  [pwB  (X,  t)  -  pwA  (X,  t)}A!b  (X)dX 
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where  the  difference  pwAB(X,t)  =  pwB(X,t )  -  pwA(X,t)  is  calculated  using  the  flow  potential 


^AB=^B-^A=~T-4(x  +  aJ)ioB[(U-aJt]  +  j-  J  log(X- a) A*b(a)d(j.  (2.34) 
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Interestingly,  the  potential  does  not  depend  on  r .  Then 
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The  wave  drag  component  associated  with  this  pressure  difference  is 
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Then  the  wave  drag  is 

D(t)  =  DA+DAB(t )  (2.36) 

Note  that  the  second  integral  in  (2.35)  is  not  singular  since  (U-aao)t<X-a<X ,  where 
(U-ax)t  <  X  <L  and  L  /( U  +  ax)<t<L  /(U  -ax).  As  t  — >  L  /(U  -ax),  the  lower  limit  behaves 
as  (U  -  ax )t  — »  L ,  and  both  integrals  tend  to  zero;  i.e.,  the  wave  drag  (2.36)  is  matched  with  DA 
at  the  upper  boundary  t  =  L/(U  -ax) . 
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2.5  Wave  drag  for  0  <t  <L  /( U  +  ax) 


In  this  case,  Region  A  corresponds  to  0  <  X  <LA  ,  Region  B  corresponds  to  LA  <  X  <  LB  ,  and 
Region  C  corresponds  to  LB  <X  <L  ,  where  LA=(U  -  ax)t  and  LB=(U  +  ax)t .  The  wave  drag 
is 
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=  Da  +  Dar  (t )  +  Dbc  (t ) 

Here  DA  is  given  by  (2.33),  DAB(t)  is  given  by  (3.35),  and  DAC(t )  is  calculated  using  the 
potential 
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The  wall  pressure  is  expressed  as 
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Finally  we  get 
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2.6  Wave  drag  for  t  — >  0 


(2.40) 


For  t  — >  0,  the  components  DM{t)  and  DBC(t )  need  a  special  treatment,  since  the  drag  tends  to 
infinity  as 


D(t)^p£L-L\4(X)fdx. 
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More  precisely,  this  limit  corresponds  to  aj/R—> qo;  i.e.,  R«aJ«L  and  pressure  waves 

propagate  at  a  distance,  which  is  much  larger  than  the  body  radius  and  much  smaller  than  the 
body  length. 

The  expression  (2.41)  can  be  derived  by  considering  (1.11)  at  e  =  aj 
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Using  (1.29)  we  write  (2.43)  as 


Then 


2  n  2  aj 


dr  2  n  r 


d<f>  U2  \  U  4(X) 

-T7  =  T-  ■ 4  W  (log  r  -  log  2a j)  -  - - b—— . 

ot  2n  2  n  t 


Using  (2.2)  we  obtain  the  wall  pressure 
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Then  the  wave  drag  is  expressed  as 
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The  first  two  terms  give 
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The  third  term  is  zero  since  Ab( 0)  =  A'h (L)  =  0  ,  and  the  forth  term  gives  the  expression 
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that  coincide  with  (2.41). 

According  to  (2.48)  the  wave  drag  tends  to  infinity  as  t  — >  0 .  To  resolve  this  singularity  we 
should  consider  a  much  shorter  time  interval  in  which  aj  =  0(R ) .  The  expression  (1.16)  gives 
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It  is  immediately  seen  that  for  very  short  time,  aj /  r  — »  0 ,  (2.49)  gives 
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Then  the  wave  drag  at  t  =  0  is 
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where  8  =  rbmsa  I L  is  the  body  thickness  ratio. 
Note  that  the  steady-flight  drag 
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Then  the  wall  pressure  is 
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and  the  wave  drag  (for  short  time  intervals)  is  expressed  as 
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The  expressions  (2.59)  and  (2.60)  are  convenient  for  numerical  integration. 

Combining  (2.59)  with  (2.40)  and  subtracting  their  common  part,  we  obtain  a  composite  formula 
for  0  <  t  <L/(U  +  ax) 
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where  the  last  term  represent  subtraction  of  the  common  part. 


2.7  Summary  of  analysis 

Summarizing  the  foregoing  results,  we  outline  formulas  for  the  wave  drag: 
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For  0  <  t  <  L/(U  +  ax) 
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These  relations  are  used  for  calculations  discussed  hereafter. 
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(2.68) 
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3.  Numerical  examples 

The  FORTRAN  code  “Wave  Drag”  was  developed  to  calculate  the  wave-drag  time  history, 
D(t) ,  using  the  relations  (2.62)-(2.68).  Numerical  examples  were  generated  for  the  body  shape 


rb{X)  =  rmsm{7tXIL). 

In  the  dimensionless  variables  X  =  X  /  L  and  7b  -rbl  L  ,  this  shape  is  expressed  as 

7b  =  8 sin(;rX)  .  (2.70) 

Calculations  were  carried  out  for  8  =  0.05 .  The  body  shape  characteristics  are  shown  in  Figure 
A4. 

D(T ) 

Figure  A5  shows  the  drag  coefficient  CD(T)  = - - — — —  for  the  case  of  Mach  number 

A^PJJ  /  2 

M  =  U / 0^=2 ,  where  T  =  taa0  / L  is  nondimensional  time  and  Amax  =  nrbmSji  =  n82l}  is 
maximum  cross-sectional  area.  The  characteristic  times  T0  =  8  (t0  =  rbwm  /  ax ),  TA  =  1/(M-1) 
( tA=L  /(U  +  ax))  and  TB  =1/(M  +  i)  (tB=  L/(U  +  ax))  are  marked  by  arrows .  The  circle  shows 
0,(0)  calculated  using  (2.69).  For  T  >TA,  the  transient  process  is  over  and  the  wave  drag 
coefficient  is  constant.  Note  that  this  coefficient  does  not  depend  on  Mach  number.  In  accord 
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with  (2.62)  it  is  a  function  of  the  body  shape  only.  The  initial  wave-drag  coefficient  CD (0)  is 
significantly  higher  than  CD  relevant  to  steady  flight.  The  drag  quickly  decreases  in  the  time 
interval  T  =  0(S)  and  then  relatively  slow  approaches  the  steady  state  level. 


Figure  A4  Characteristics  of  the  body  shape  (2.70),  5  =  0.05 . 


Figure  A5  Wave  drag  coefficient  CD(T )  for  M  =  2 ,  5  =  0.05 . 
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Figure  A6  shows  CD(T)  for  different  Mach  numbers.  As  expected,  the  transient  time  increases 

and  tends  to  infinity  as  the  Mach  number  decreases  and  approaches  M  =  1  (transonic  case).  Note 
that  the  transonic  case  needs  special  treatments,  since  the  linear  acoustic  theory  is  not  valid  in  the 
vicinity  of  M  =  1  even  for  slender  bodies. 


Figure  A6  Wave  drag  coefficients  CD(T )  for  various  Mach  numbers,  5  =  0.05 . 


4.  Comparison  with  Euler  CFD  solution 

Numerical  calculations  of  an  axisymmetric  inviscid  flow  over  a  slender  body  of  revolution  were 
carried  out  for  the  free-stream  parameters: 

•  Mach  number  Afx  =  2 

•  Static  pressure  px  =  383410  Pa 

•  Temperature  Tx  =  277  K 

•  Density  px  =4.8222  kg/m3 

•  Velocity  =  667.03  m/s 

The  body  shape  is  given  by  (2.70)  with  the  parameters: 

•  Body  length  L  =  1  m 

•  Thickness  ratio  5  =  0.01 
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Near  the  nosetip  X  — >  0 ,  the  body  radius  is  rb  ~  8k X  that  gives  the  nose  half-angle 

a  »  8k  =  =  1.8°.  Note  that  the  Mach  wave  angle  is  y  =  arcsin  ( M”1 )  =  30°  . 

100  6  v  ' 


T=ta/L 


Figure  A7  The  drag  coefficient  CD  as  a  function  of  time  T  =  taao/  L ,  solid  line  -  CFD,  symbols 
-  theoretical  solution  summarized  in  Section  2.7. 


The  computational  grid  contains  approximately  105  cells:  «300  grid  points  along  the  body 
surface  and  «  300  grid  points  in  the  direction  normal  to  the  body  surface. 

The  drag  coefficient  versus  time  is  shown  in  Figures  A7  and  A8  (black  lines)  for  relatively  large 
and  short  initial  time  intervals,  respectively.  The  wave  drag  almost  instantaneously  reaches  its 
maximum,  decreases  with  time  to  the  minimum  CD  ®  3.65  xlO-3  at  T  «  0.13  ,  and  then  tends  to 

the  steady-state  level  CD  « 4.4x1 0‘3 .  The  CFD  solution  (black  line)  agrees  well  with  the 
theoretical  solution  summarized  in  Section  2.7. 

Instantaneous  pressure  fields  at  various  time  moments  (CFD  solution)  are  shown  in  Figure  A9. 
In  accord  with  linear  acoustic  theory,  the  body-induced  pressure  waves  propagate  along 
characteristics  and,  finally,  form  a  steady-state  pressure  field  around  the  body  moving  with 
constant  Mach  number.  For  this  steady  field,  distributions  of  axial  and  radial  velocity, 
temperature  and  Mach  number  are  shown  in  Figures  A10-A13. 
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T=ta/L 


Figure  A8  The  drag  coefficient  CD  as  a  function  of  T  =taao/  L  for  the  initial  time  interval,  black 
line  -  CFD,  red  line  -  theoretical  solution  summarized  in  Section  2.7,  red  circle  - 
Co(0)  calculated  using  (2.69). 
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Figure  A9  Evolution  of  static  pressure  field. 
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Figure  A10  Axial  velocity  at  t  =  3000  x  10’6  sec. 
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Figure  All  Radial  velocity  at  t  =  3000  x  10'6  sec. 
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Figure  A12  Temperature  at  t-  3000  x  10‘6  sec. 
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Figure  A13  Mach  number  at  t  =  3000  x  10"6  sec. 
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5.  Conclusions 

1 .  Mathematical  modeling  of  unsteady  effects  associated  with  acceleration  of  the  body  center  of 
gravity  has  been  carried  out.  Namely,  we  have  conducted  theoretical  analysis  of  the  wave 
drag  generated  by  a  non-lifting  body  of  revolution  during  its  instantaneous  supersonic  start. 
The  analysis  is  based  on  linear  acoustic  theory  and  asymptotic  techniques. 

2.  Typical  phases  of  the  transient  process  were  identified.  Analytical  expressions  of  the  flow 
potential  were  derived  for  each  phase.  These  expressions  were  used  to  obtained  the  pressure 
distribution  on  the  body  surface  and  calculate  the  wave  drag. 

3.  It  was  shown  that  in  the  initial  time  interval  R « aj « L ,  the  wave-drag  solution  is 

singular.  Asymptotic  analysis  of  this  initial  phase  of  transient  process  allowed  us  to  resolve 
the  singularity,  obtain  an  analytical  expression  for  the  initial  wave  drag  D(t  =  0)  and  build 
up  a  composite  solution  convenient  for  computations. 

4.  The  aforementioned  theoretical  results  were  incorporated  into  a  FORTRAN  code  providing 
quick  calculations  of  the  wave  drag  time  history.  Numerical  examples  were  generated  for  a 
slender  body  of  half-sine  shape  with  thickness  ratio  5  =  0.05  for  different  Mach  numbers 
(M=1.2,  1.4,  1.6  and  2)  of  instantaneous  body  start.  It  was  shown  that  the  initial  wave-drag 
coefficient  CD(t  =  0)  is  significantly  larger  than  CD  relevant  to  steady  flight.  The  wave  drag 

quickly  decreases  in  the  time  interval  t  ~  SL!  ax  and  then  relatively  slow  approaches  the 

steady-state  level.  The  transient  time  increases  and  tends  to  infinity  as  the  Mach  number 
decreases  and  approaches  M  =  1  of  transonic  flow. 

5.  The  transient  process  was  also  simulated  numerically  by  solving  3-D  Euler  equations  for 
Mach=2  start  of  the  aforementioned  body  with  S  =  0.01 .  It  was  shown  that  the  theoretical 
wave  drag  D(t)  is  in  excellent  agreement  with  that  predicted  by  CFD  in  all  phases  of  the 
transient  process. 

These  results  provide  a  good  launching  pad  for  theoretical  modeling  of  unsteady  body  motions 
including  the  coupling  between  body  dynamics  and  aerodynamics. 
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